Neutron scattering study of crystal fields in CeRhIn 5 
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Neutron scattering results for the tetragonal compound CeRhlns give evidence for two crystal field 
(CF) excitations at 6.9 and 23.6 meV. The scattering can be fit assuming a set of CF parameters 
B2 = -1.03 meV, B4 = 0.044 meV and B 4 = 0.122 meV. To compare our results to previous work, 
we calculate the susceptibility and specific heat for this CF scheme, including a molecular field term 
A =35 mol/emu to account for the Kondo effect. We also include a calculation based on these CF 
parameters that uses the non-crossing approximation to the Anderson model to estimate the effect 
of Kondo physics on the susceptibility, specific heat and neutron linewidths. 

PACS numbers: 75.30.Mb 75.20.Hr 71.27.+a 71.28.+d 61.10.Ht 



CeRhlns crystallizes in the same tetragonal HoCoGas 
structure as the_heavy fermion superconductors Celrlns 
and CeCoIngtHj'H. At ambient pressure CeRhlns under- 
goes a transition to an antiferromagnetic (AF) state at 
IV = 3.8 KBq. With application of hydrostatic pres- 
sure the Neel temperature remains essentially constant 
until antiferromagnetism disappears and superconduc- 
tivity appears _at pressures above 15 kbaiu. Recently, 
Pagliuso et alB have suggested the importance of CF 
splitting to the ground state properties of the CeMIn 5 
family of heavy fermion superconductors, underscoring 
the fact that the ultimate ground state achieved by a 
particular member of the family must grow out of the 
ground state crystal field doublet. Thus a careful deter- 
mination of both the CF splitting and wavefunctions is 
important. To that end we have begun to directly probe 
the CF energy level splitting in the CeMIns family using 
inelastic (IE) neutron scattering. The first step in our in- 
vestigations has been determining the crystal field level 
scheme in CeRhlns. 

In CeRhlns, as in the other members of CeMIns family, 
the crystal field Hamiltonian in tetragonal symmetry can 
be written 



H, 



CF 



S 9 O 9 ° + B Q A O A + BtOi 



where OJ™ and B\ 



are the Stevens operators and CF pa- 
rameters respectively. The Ce 3+ J = 5/2 wavefunction 



splits into three doublets, T 



(i) 



{a[±5/2> + [3[i 



{/3[±5/2) - a[±3/2)} and T e = [±1/2)1: 



: 3/2)}, 



analysis of susceptibility and thermal expansion resultscl 
suggested crystal field levels Tj 1 ^ and Tq at £'=0, 

5.86 meV (68 K) and 28.43 meV (300 K) respectively, 
with j3 — 0.969 (yielding a nearly pure [±5/2) ground 
state). A subsequent studyO based on an analysis of 
the susceptibility and specific heat suggested a similar 
scheme, but with splittings 6 and 12 meV (70 and 140 
K). In this paper we report the results of an analysis of 



neutron scattering data for CeRhlns which indicate that 
these initial estimates are nearly correct; our results have 
somewhat different values for the splittings and a smaller 
value for the mixing parameter (3, i.e., a greater admix- 
ture of [=p3/2) into the [±5/2) ground state. To assi: 
in comparison of our results to those of Pagliuso et 
and Takeuchi et alu, we report calculations of the specific 
heat and magnetic susceptibility based on our CF param- 
eters which include the Kondo effect in an ad hoc manner 
similar to those of refs. 5 and 6. We also present more 
sophisticated calculations that employ the non-crossing 
approximation (NCA)LrtJ to the Anderson model in order 
to estimate the effect of Kondo spin fluctuations on the 
susceptibility, specific heat and IE neutron spectra. 

Large high quality single crystals of CeRhlns and 
LaRhlns were obtained using the flux-growth methodtl. 
For the magnetic susceptibility and specific heat mea- 
surements, single crystals were carefully prepared which 
were free of residual In flux; in the case of the neutron 
scattering measurements, ~50 g of single crystals for both 
CeRhlns and LaRhlns were powdered. The neutron scat- 
tering experiments were performed in time-of-flight mode 
using LRMECS at IPNS (Argonne National Laboratory) 
with experimental conditions that were similar to that 
of an earlier reports!. A key problem in our investiga- 
tions was the high neutron absorption of both In and 
Rh. In initial experiments the standard LRMECS sam- 
ple holder was used; however, in subsequent experiments 
a new sample holder was employed which was designed to 
maintain a more uniform sample thickness than the stan- 
dard holder, thus allowing for a more accurate absorption 
correction. Neutron scattering spectra were collected for 
several different incident energies (Ei) and temperatures 
between 8 and 140 K with counting times ranging from 
24 to 48 hours. To improve statistics, we were able to 
take advantage of the nondispersive nature of the CF 
scattering and group detectors into three bins with mean 
scattering angle 20° (low Q), 60° and 100° (high Q). A 
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FIG. 1: Neutron energy spectra of (a) CeRhlns and (b) 
LaRhlns at an initial energy Ei = 35 meV, at 8 K and for two 
mean scattering angles, 20° and 100°. The data have been 
corrected for neutron absorption and the scattering from the 
sample holder has been subtracted from the data, (c) The 
Q=0 magnetic scattering, determined as described in the text, 
in CeRhlns at 8K and for three incident energies Ei. 



Vanadium standard was utilized to put the scattering on 
an absolute scale. 

Data for CeRhlns and LaRhlns (measured to help 
identify the nonmagnetic scattering in CeRhlns) at 8 K 
and Ei = 35 meV for low and high Q are shown in Fig. 
[]]. The data were corrected for absorption assuming a 
uniformly thick flat-plate sample and the scattering due 
to the empty sample holder was then subtracted from the 
data. Direct comparison of low angle scattering (where 
magnetic scattering is strongest) for CeRhlns (Fig. 0a) 
and LaRhlns (Fig. |l|b) shows two additional peaks near 
7 and 23 meV. In particular, we determine the nonmag- 
netic scattering in CeRhlns in two ways: 1) By using 



10 



_ 

> 

10 



- 5 

c? 



CO 



10 







(a) 



T~ji 1 1 1 1 1 1 1 ' 

T = 8 K 

CeRhln 

5 



T = 70K 




= 140 K 








10 20 

AE (meV) 



30 



FIG. 2: Temperature dependence of the magnetic part of the 
IE neutron scattering response of CeRhlns for Ei — 35 meV. 
The scattering dependence due to the Ce 3+ form factor has 
been removed as in the previous figure. The data at all three 
temperatures ((a) 8 K, (b) 70 K, and (c) 140 K) have been fit 
simultaneously (solid lines) with a least squares fitting routine 
to determine the crystal field parameters. The results of the 
fitting parameters including the crystal field parameters are 
displayed in table We have included in (a) the results of 
the NCA calculation (dashed line). 



the expression S mag (20°) = S(Ce,20°) - fS(La,20°) 
where we choose the factor / as the ratio (0.75) of the 
total scattering cross-sections cr(CeRhIns)/cr(LaRhIn5). 
2) By determining the ratio R = S(La, 100°) /S (La, 20°) 
for scaling the high angle nonmagnetic scattering to low 
angletL3. Excellent agreement with 1) is obtained using 
S mag (20°) = S{Ce, 20°) - FS(Ce, 100°)/R with inclu- 
sion of an additional factor F = 1.33 to account for the 
difference in Q-scaling of the La and Ce compounds. The 
value of E is similar to the one used in a recent study of 
YbXCuJj; it can be justified on the basis that for high an- 
gle scattering the data are predominantly single-phonon, 
proportional to a, while the low angle scattering contains 
a significant contribution from multiple scattering (one 
elastic and one phonon) proportional to a 2 , so that the 
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TABLE I: Crystal field parameters B™ , splittings and 
Lorentzian halfwidths F of the IE excitations at four temper- 
atures for CeRhIn5 and the wave function mixing parameter 
j3. The units of all quantities (except for /3, which is unitless) 
are meV. The reduced Chi-square for the fit was \ 2 = 0-69 
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cross section does not cancel in the ratio. Results of this 
analysis for three different Ei are shown in Fig. lc. The 
dependence of the scattering on the Ce 3+ form factor has 
been removed in this plot, so the data represent the Q = 
scattering with the assumption that the crystal fields 
are in fact purely local and uncoupled entities. The data 
have been truncated below 0.15-Ei (where the elastic line 
dominates the scattering) and above 0.8£i, where statis- 
tics are small due to the kf/ki factor. Good agreement is 
evident for data taken at three different Ei, with all data 
sets displaying magnetic excitations at approximately 7 
and 24 meV. 

In Fig. 2 we plot the Q = (form factor removed) mag- 
netic scattering (method 1), determined at Ei — 35 meV, 
for three different temperatures. We have performed a si- 
multaneous least squares fit to four datasets (8 K, 70 K 
and 140 K at Ei = 35 meV and 8 K at Ei = 80 meV) 
to determine the CF parameters. The fit includes the 
effects of instrumental resolution. Variables of the fit in- 
clude i?", B®, B\ and an overall scale factor (which four 
parameters were constrained to the same values for all 
datasets) and the Lorentzian halfwidth T of the IE ex- 
citations which was allowed to vary with temperature. 
(We constrained the quasi-elastic (QE) halfwidth to 1/2 
F.) Results of the fit are shown in Table | and plotted in 
Fig. 2. 

To compare our results to those of Pagliuso et al.a and 
Takeuchi et alu, we have calculated the susceptibility 
and specific heat (Fig. 3). The susceptibility includes 
a positive molecular field contribution A = 35 mol/emu 
where A represents contributions to 1/x from AF and 
Kondo fluctuations. At high T these contribute to 1/x as 
(T K + T N )/Cj; with Cj = 0.807 emu-K/mol for J = 5/2 
and T N = 3.8 K this gives T K ~ 25 K. We note that this 
value of UbTk is similar to the width of the 7 meV IE 
excitation at 8 K. The calculation for the specific heat 
contains both a Schottky te«n due to the excited levels 
and a Kondo doublet termlill with Tk — 25 K for the 
ground state level, which puts the calculated specific heat 
in the range 20-50 K in better agreement with measured 
value - without this, the calculated value due only to the 
Schottky contribution is smaller by a factor of 0.8. We 



have not attempted to fit for the effects on C mag and \ 
of the AF transition at 3.8 K. 
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FIG. 3: (a) Measured anisotropic susceptibilities \ zz (tri- 
angles) and \ xx (circles) for CeRhlns compared to the value 
calculated for the CF parameters of Table @ with a molecu- 
lar field contribution A = 35 mol/emu (solid lines) and com- 
pared to the results of the NCA calculation (dashed lines), (b) 
Magnetic specific heat compared to the value calculated for 
a Schottky contribution from the excited levels and a Kondo 
contribution from the ground state doublet (solid line) and to 
the results of the NCA calculation (dashed line). 

A more sophisticated way to include Kondo spin fluctu- 
ations is through calculation for the Anderson model. We 
present resultsjabtained using the non-crossing approxi- 
mation (NCAjLra. We have used a Gaussian background 
band with density of states N(e) = e^ e / w ^ /^wW with 
W = 3 eV and we set the 4/ level position at Et = 2 eV 
and the spin-orbit splitting at E so — 0.273 eV, which 
are standard values for Ce. Since the Kondo physics 
renormalizes the CF levels upward by an amount approx- 
imately equal to the Kondo temperature the bare level 
energies were chosen to be Eb = 5.3 meV and E c = 23 
meV which are smaller than the measured level energies. 
The mixing parameter [3 — 0.80 was chosen to be similar 
to that obtained in Table 1. The hybridization was then 
varied until a good fit to the anisotropic susceptibility 
was obtained for V — 0.4665 eV. The results for S mag , \ 
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and C mag are given in Figs. 2a and 3. 

We now turn to discussion of the effect of system- 
atic errors on our conclusions. As mentioned previously, 
the neutron absorption of In and Rh is an important 
consideration. Comparison of the data for two differ- 
ent sample holders (which exhibited small differences in 
sample thickness and distribution) indicated similar re- 
sults, augmenting our belief that the absorbtion correc- 
tion employed is correct. If the nonmagnetic background 
subtraction is varied by varying / (method 1) or F/R 
(method 2), the scattering at 7 meV is relatively unaf- 
fected but the strength of the 24 meV scattering, and 
hence (3, is affected somewhat. Given the good consis- 
tency between results at different Ei and T we think that 
our CF scheme is basically correct. We were unable to 
observe quasi-elastic (QE) scattering, due to the require- 
ment that to obtain the resolution necessary a small Ej 
is required which causes the effects of neutron absorp- 
tion, which varies as 1/y/E, to become large. In our fits 
we constrained the QE halfwidth to half the value of the 
IE width to prevent proliferation of fit parameters. Con- 
straining to other values (e.g. Tqe = Tje) leads only to 
minor variation in the final fits. 

Our fits to x using the CF parameters plus a molec- 
ular field term are nevt as good as those of Takeuchi et 
al.U or Pagliuso et al.U. However, their fits use a value 
of (3 very close to unity, indicating essentially no [=p3/2) 
admixture into [±5/2) ground state. In this case there 
would be no observable amplitude for the Am 2 = 1 tran- 
sition to [±1/2) state at 24 meV. This cannot be correct 
as we clearly observe this transition in the neutron scat- 
tering data. A possible reason that our fits are not as 
good as those of refs. 5 and 6 is that we do not include 
the effect of exchange anisotropy, which should only be 
important below 20.-K. Such anisotropy can be mimicked 
as in Pagliuso et alB through inclusion of an anisotropic 
mean field parameter, which we have chosen not to do 



for simplicity. 

On the other hand, the NCA calculations based on our 
CF scheme and a Kondo temperature of order 25 K does 
an excellent job reproducing x- However, it overestimates 
the width of the 7 meV excitation as seen in Fig. 2a and 
underestimates the temperature of the peak in the spe- 
cific heat (Fig. 3b) . These deviations from the data may 
reflect the fact that we have neither included antiferro- 
magnetic exchange, exchange anisotropy nor anisotropic 
hybridization (i.e. different hybridization to the different 
CF multiplets) in the NCA fits. 

In summary, we find a more significant [=p3/2) admix- 
ture into [±5/2) (/3 = 0.80) ground state than found 
earlier by Takeuchi et al. (/3 = .969)0 or Pagliuso et 
al. ((3 ~ 1)0. The resulting CF level parameters pro- 
vide reasonable fits to both the magnetic susceptibility 
and specific heat with the inclusion of a mean field pa- 
rameter and a Kondo doublet respectively. In addition, 
NCA fits the susceptibility remarkably well with some 
deficiencies in both the specific heat and neutron scat- 
tering linewidths. Taken together the NCA calculations 
and the fits to specific heat and susceptibility all indi- 
cate a Tk ~ 25 K. We note that the ordered moment 
g^B(Jz) = 0.92/zb deduced for [3 =0.80 is substantially 
larger than the value 0.36/j,n needed to fit the diffraction 
pattern in the ordered statea. For our estimate of Tk the 
moment is reduced by the Kondo physics at temperatures 
T>T N . 
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